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Abstract 



> 
o 

I The BaBar Collaboration has recently reported the measurement of the ratio of the 

branching fractions of B — > D(D*)t~v t to B — > D(D*)£~v e which deviates from 
the Standard Model prediction by 2a(2.7cr). This deviation goes up to 3. 4a level 
■ when the two measurements in the D and D* modes are taken together and could 

indicate new physics. Using an effective Lagrangian for the new physics, we study 
the implication of these results and calculate other observables that can shed light 
^ ; on the nature of the new physics. We show that the measurements of the forward- 

backward asymmetries and the r and D* polarization fractions can be distinguished 
among the various couplings of the new physics operators. 

1 Introduction 

The Standard Model (SM) has been extremely successful in furthering the under- 
standing of the various measurements of Branching Ratios (BR) and asymmetries 
in the quark sector. In the quark flavor sector, the B factories, BABAR and Belle, 
have produced an enormous quantity of data in the last decade. There is still a lot of 
data to be analyzed from both experiments. The B factories have firmly established 
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the CKM mechanism as the leading order contributor to CP violating phenomena in 
the flavor sector involving quarks. New physics (NP) effects can add to the leading 
order term producing deviations from the Standard Model (SM) predictions. But 
even after this incredible success the mechanism of Electroweak Symmetry Breaking 
(EWSB) and the origin of masses in the SM still remain very poorly understood. 
In this respect, the second and third generation quarks and leptons are quite spe- 
cial because they are comparatively heavier and are expected to be relatively more 
sensitive to new physics. As an example, in certain versions of the two Higgs dou- 
blet models (2HDM) the couplings of the new Higgs bosons are proportional to the 
masses and so new physics effects are more pronounced for the heavier generations. 
Moreover, the constraints on new physics involving, specially the third generation 
leptons and quarks, are somewhat weaker allowing for larger new physics effects. 
Interestingly, the branching ratio of B — > tv t shows some tension with the SM pre- 
dictions m and this could indicate NP [2], possibly coming from an extended scalar 
or gauge sector. There is also a seeming violation of universality in the tau lepton 
coupling to the W suggested by the Lep II data which could indicate new physics 
associated with the third generation lepton [3]. 

If there is NP involving the third generation leptons one can search for it in B 
decays such as B — > Dt~i7 t1 B — > D*t~i7 t [4]. The semileptonic decays of B meson 
to the r lepton is mediated by a W boson in the SM and it is quite well understood 
theoretically. In many models of new physics this decay gets contributions from 
additional states like new vector bosons or new scalar particles. The exclusive decays 
B — > Dt~i7 t and B — > D*t~i7 t are important places to look for NP because, being 
three body decays, they offer a host of observables in the angular distributions of 
the final state particles. The theoretical uncertainties of the SM predictions have 
gone down significantly in recent years because of the developments in heavy-quark 
effective theory (HQET). The experimental situation has also improved a lot since 
the first observation of the decay B — > D*r~v T in 2007 by the Belle Collaboration [5]. 
After 2007 many improved measurements have been reported by both the BaBaR 
and Belle collaborations and the evidence for the decay B — > Dt~i7 t has also been 
found [HI El IB] • Recently, the BaBar collaboration with their full data sample of an 
integrated luminosity 426 fb _1 has reported the measurements of the quantities [H] 

0.440 ±0.058 ±0.042, 

= 0.332 ±0.024 ±0.018. (1) 

The SM predictions for R(D) and R(D*) are QUI EI] 

R(D) = 0.297 ±0.017, 

R(D*) = 0.252 ±0.003, (2) 



R{D) 
R(D*) 



BR(B ->■ Dt~u t ) 
BR(B -> D£-u e ) 
BR(B ->■ D*t-u t ) 
BR(B -> D*i-v t ) 
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which deviate from the BaBar measurements by 2a and 2.7a respectively. The 
BaBar collaboration themselves reported a 3Aa deviation from SM when the two 
measurements of Eq. [1] are taken together. 

These deviations could be sign of new physics and already certain models of new 
physics have been considered to explain the data [T2]. In this work, we calculate 
various observables in B — > Dt~v t and B — > D*t~i7 t decays with new physics. We 
write the most general effective Lagrangian that affect these decays. The Lagrangian 
contains two quark and two lepton scalar, pseudoscalar, vector, axial vector and 
tensor operators. Considering subsets of the NP operators at a time, the coefficient 
of these operators can be fixed from the BaBar measurements and then one can 
study the effect of these operators on the various observables. 

The paper is organized in the following manner. In Sec. 2 we set up our for- 
malism where we introduce the effective Lagrangian for new physics and define the 
various observables in B — > Dt~v t and B — > D*t~i7 t decays. We also present the 
SM predictions for these observables in that section. In Sec. 3 we present the nu- 
merical predictions which include constraints on the new physics couplings as well 
as predictions for the various observables with new physics. Finally, in Sec. 4 we 
summarize the results of our analysis. 

2 Formalism 

In the presence of new physics (NP), the effective Hamiltonian for the quark- level 
transition b — > cl~vi can be written in the form [13] 



^^±[(1 + V L ) [c^P L b] [h"P L vi\ +V R [crfP R b] [h^Pm] 
+S L [cP L b][lP L u l ] +S R [cP R b] [IPm] +T L [ca^P L b] [Ta^PLvA (3) 



where Gp = 1.116637 x 10~ 5 GeV~ 2 is the Fermi coupling constant, V c b is the 
Cabibbo-Koboyashi-Maskawa (CKM) matrix element. Pl,r = (1 TTs)/^ is the 
projector of negative /positive chirality, and we use <t m „ = i[ r y fJ ,,j v ]/2. We have as- 
sumed the neutrinos to be always left chiral. Further, we do not assume any relation 
between b — > ul~vi and b — > cl~vi transitions and hence do not include constraints 
from B —7- tv t . The SM effective Hamiltonian corresponds to Vl = Vr = Sl = 
Sr = Tl = Tr = 0. In this paper we will ignore the tensor interactions. With this 
simplification we write the effective Lagrangian as 

G V 

Heff = F ^ b { r^ 1 ~ 75 ) & + dvcif.b + ^7^75^7^(1 - Ts)^ 

+ g s cb + gpc^b 1(1 - 75)^ + h.c}, (4) 
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where g v ,A = V R ±V L and g s ,p = S R ± S L . 
We will now consider the two cases: 



Case a : In this case we will set Sl, Sr = and assume that the NP affects lep- 
tons of only the third generation. This scenario could arise from the exchange 
of a new charged W boson [TJ 



Case b : In this case we will set Vl,Vr = and assume that the NP affects 
only leptons of the third generation. This scenario could arise in models with 
extended scalar sector [151 . 



The polar angle differential decay distribution in the momentum transfer squared 
q 2 for the process B — > D^Hui can be written in the form 

Tir^T = E \M{B ^ D^l n )\ 2 , (5) 

dq 2 d cos 9 l 2567r 3 m| pola ^ ation 



where V\ = y 1 — mf /q 2 and the momentum of the meson in the B meson rest 
frame is denoted as \p D (*)\ = \ l ^ 2 {m 2 B , m 2 DM , q 2 )/2rriB with A(a, b, c) = a 2 + b 2 + c 2 — 
2(ab + be + ca). 

2.1 B — > D*r~i7 T angular distribution 

The full B — > D*t~i7 t angular distribution is given by, 

dY D * 



dq 2 d cos 9i 



N\p D .\[2\Ao\ 2 sin 2 e l + {\A\\\ 2 + |^ ± | 2 )(1 + cos^) -^e^n^l] cos^ 

2 

+ ^ (2\Ao cos 6 l - A tP \ 2 + (|^|| | 2 + \A±\ 2 ) sin 2 ei)] , (6) 



where Q[ is the angle between the D* meson and the r lepton three-momenta in the 
q 2 rest frame, N = (l — and the amplitude A t p is 

Ap = Ut + ^-Ap). (7) 

The differential decay rates for the r helicities, A r = ±1/2, various transversity 
amplitudes and the form factors are defined in appendix [B] and [D] respectively. 

The angular distribution allows us to define several observables pUl ELL]. The 
starting point is to obtain the decay rates dT/dq 2 for the r helicities, A T = ±1/2 , 
after performing integration over cos 9i , 
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dr D *[\ T = -1/2] 

dq 2 

dr D *[\ T = 1/2] 

dq 2 

where |Ar| 2 = |A| 2 + |*<4||| 2 + \A± 
differential branching ratio (DBR): 



8N\p 



D\ 



3 

4N\p D \ 



q 2 



Ut| 2 + 3|Ap| 2 



The summation of these rates give the total 



dBr[B D*t-i7 t ] 
dq 2 



SN\ Pd \t b 

3h 



m 
2q 2 



3m; 
~2f 



tp\ 



(9) 



where tb is the B meson life-time. Furthermore, one can also explore the q 2 depen- 
dent ratio 



dBr[B D*T-i7 T ]/dq 2 
dBr[B DH-v^jdq 2 



(10) 



where I denotes the light lepton (e, //). The ratio Rr,*(q 2 ) is independent of the form 
factor h^iw). The SM predictions of DBR and Ri>*(q 2 ) for the decay B — > D°*t~v t 
are shown in Fig. [TJ The numerical values of B — > D* form factor parameters 
^1(1)1^1, p 2 and #0,1,2(1) are given in appendix ID1 In the SM, the DBR for the 
decay B D°*r~u T peaks (« 0.27% GeV" 2 ) at q 2 « 8 GeV 2 , and the ratio # D *(g 2 ) 
rises to more than 50% at large q 2 . It is clear from the plot that the uncertainty in 
Rd*{q 2 ) is less than that in the DBR. 
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0.00 
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Figure 1: The left (right) panels of the figure show the q 2 dependence of DBR 
(i?D*(g 2 ) ) for the decay B — > D°*t~u t . The bands correspond to uncertainties in 
^4i(l)|Kfc|, P 2 and #0,1,2(1)- The errors are added in quadrature. 

Next, we define the forward-backward asymmetry (AFB) in the angular distri- 
bution by integrating over cos Q\ as 



[A FB }D*(q 2 



(Jo ~ J°!)d cos 9i 



dV 



D* 



dq 2 d cos di 



dT D * 
dq 2 



3 (RejAuA l] + ^Re[A A, 
2 \A T \ 2 (l 



tp. 



■2c,- 



if 



The perpendicular transversity amplitude A± is proportional to \/w 2 — 1 (see ap- 
pendix |B] for the details), hence for the light leptons [A FB ] D *(q 2 ) vanishes at the 
end-points due to the kinematics. One can obtain the angular distribution only 
for the transversely polarized D* meson from Eq. ([6]) by dropping the amplitudes 
Ao and Atp- We now define the forward-backward asymmetry for the transversely 
polarized D* meson by integrating over cos#; as in Eq. ( jTTj) [16J: 

3 Re[A\\A* A 

~^{\A\\\ 2 + \A i 



[A T FB W(q 2 ) 



(12) 



Fig. [2] shows the SM predictions for [A FB ] D *(q 2 ) and [A^ B ] D *{q 2 ) in the decay 



B ->■ D*r 

J2 



[A FB \ D *(q 2 

t2 



is ~ 20% and negative at low q 2 , and has a zero crossing 
at q z ~ 5.64 GeV^. However, the asymmetry [Aj^ B ]D*(q 2 ) is always positive and 
large (~ 40%) at low q 2 . 





Figure 2: The left (right) panel shows the q 2 dependence of [Ap#]d* ([^4^ b ]d* ) for 
the decay B — > D *t~v t . The bands correspond to uncertainties in p 2 and _R ,i,2(l)- 
The errors are added in quadrature. 

We now define the longitudinal and transverse polarization fractions of the D* 
meson from Eq. M as 



FT 



IA| 2 (i + ^) + g|Apf 
\A T \*(l + ?$) + ^\A tP \ 



(\A\\\ 2 + \A ± 


I 2 )(l + 


rr&\ 
2? J 




\A T 


2 ( 1 + #) 


, 3m 2 


A tP 


12 



(13) 



where F®* + F F * = 1. The D* polarization fractions can be measured by fitting to 
the decay distribution in Eq. [6] or from D* decays. 

Finally, one can also define the longitudinal polarization fraction of the r lepton 
in the q 2 rest frame as 



dT D *[\ T = 


-1/2] dr D *[X T =l/2] 


\A T \ 2 (l- 


m 2 \ 


3m? 


\A tP \ 2 


dq 2 


dq 2 


2q 2 ) 


2q 2 




dT D * 
dq 2 


' \A T \ 2 (l + 


m 2 \ 
2q 2 ) 


, 3m? 
2q 2 


A tP \ 2 



(14) 



The r polarization can be measured from the decays of the r. 

The polarization fractions are independent of the form factor hA 1 (w). The SM 
predictions of the longitudinal polarization fractions of D* and r are shown in Fig. [3j 
In the SM, (q 2 ) can be as large as 0.75 at low q 2 , and it decreases to about 0.4 
at high q 2 . On the other hand, Pl T (q 2 ) is about -0.18 at very low q 2 and has a 
zero-crossing at q 2 « 3.64 GeV 2 and increases to about 0.7 at high q 2 . 




« 2 [GeV 2 ] 9 2 [GeV 2 ] 

Figure 3: The left (right) panel shows the q 2 dependence of F^* (P£ T ) for the decay 
B — > D°*t~v t . The bands correspond to uncertainties in p 2 and i?o,i,2(l)- The 
errors are added in quadrature. 



2.2 B — > Dr~v r angular distribution 

The full B — > Dt~i7 t angular distribution can be written as, 

dT D r m 2 i 

N\p D \ 2\H Q \ 2 siv 2 6 l + 2^ r (H oC os9 l -H tS ) 2 , (15) 



dq 2 d cos 9i L g2 

where 



Htg = (Ht-^-Hs). (16) 



The differential decay rates for the r helicities, A T = ±1/2, and helicity amplitudes 
H are defined in appendix O 

As in the previous section, we can define several observables using the B — > 
Dt~v t angular distribution pj]. The starting point is to obtain the decay rates 
dT/dq 2 for the r helicities, A T = ±1/2. After performing integration over cos 6i, one 
can obtain: 



dT D [X T = -1/2] 



dq 2 3 



N\p D \\H 



2 



Dl\ _ 1 /ol A m 2 



dT [X T = 1/2 



= -iV|^|-^(|if |^ + 3|^5| 2 ), (17) 
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and the summation of these differential decay rates give the DBR 



dBr[B -> Dt-j7 t ] 



SN\ Pd \t b 



ts\ 



As in the previous section, one can also explore the q 2 dependent ratio 

2 dBr[B -» Dr-j7 T }/dq 2 

^ D(g) = dBr[5 D£~Ui\/dq 2 ' (19) 

The ratio Ro{q 2 ) is independent of the form factors. The SM predictions of DBR 
and Ro{q 2 ) for the decay i? — > D°t~v t are shown in Fig. HI The numerical values 
of the free parameters in the B — > D form factors are given in appendix O In the 
SM, the DBR for the decay B ->■ D q t~V t peaks (« 0.14% GeV~ 2 ) at q 2 « 7 GeV 2 , 
while Rr,(q 2 ) shows almost a linear behavior with g 2 . 




6 8 10 

? 2 [GeV 2 ] 



Figure 4: The left (right) panel shows the q 2 dependence of DBR ( R£>{q 2 ) ) for the 
decay B — > D°t~v t . The bands correspond to uncertainties in Vi(l)|T4&|> P 2 and p\. 
The errors are added in quadrature. 

Next, we define the forward-backward asymmetry in the angular distribution by 
integrating over cos Q\ as 



[A FB \D{q 2 



U-i -Jo)dcos6, 



dq 2 d cos t 



3mt 



Re[H H t * s ] 



dT D 
dq 2 



2q 2 



\H \ 2 (l + " 



2 , 
2? 



i 3mg | rr 12 



(20) 



Finally, we define the longitudinal polarization fraction of r in the q 2 rest frame as 



dr°[A T =i/2] 

dq 2 



dr D [\ T =-i/2] 

dq 2 



, 3m£ I tt I 



dT D 
dq 2 



\Ho\ 



'1 _l_ us.' 

^ + 2q 2 ' 



3m 2 | rr |2 

-^I-Htsl 



(21) 



The r polarization can be measured from the decays of the r. 

Fig. [5] shows the SM predictions for [Ap^] £>(g 2 ) and P[(q 2 ) for the decay B — » 
D°t~v t . The forward-backward asymmetry, [A x 7 B ] Z) (g 2 ), is about ~ 50% at low q 2 
and decreases with increasing q 2 . The r polarization, P£, is about 0.4 at low q 2 and 
starts to increase after q 2 = 8GeV 2 . 

In the next section we shall study the effect of the new physics couplings Vl,r, 
and Sl r on the above observables. 
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Figure 5: The left (right) panel shows the q 2 dependence of [Apb]d (P[) f° r the 
decay B — > D°t~v t . 

3 Numerical analysis with NP 

In the numerical analysis, as indicated earlier, we consider two cases to study the new 
physics effects on DBR, the ratios Ro(d*){<1 2 ), the forward-backward asymmetries, 
and polarization fractions. In the first case we consider only vector/axial- vector 
NP couplings while in the second case we consider only scalar /pseudoscalar NP 
couplings. The numerical values of B — > D and B — > D* form factors in the heavy 
quark effective theory framework are summarized in appendix [D] A detail discussion 
of these form factors can be found in [T7] . 

In our numerical analysis, we constrain both complex/real NP couplings Vl,r 
and Sl,r using the measured R(D) and R(D*) in Eq. (JTJ at 95% C.L. We also vary 
the free parameters in the form factors discussed in appendix |D] within their error 
bars. All the other numerical values are taken from [18] and [19] . The allowed ranges 
for NP couplings are then used for predicting the allowed ranges for the observables 
discussed earlier. 

3.1 B -> D°*t-is t 

3.1.1 Pure Vl and Vr couplings present 

The combination of the couplings gy = Vr + Vl appears in both R{D) and R(D*), 
while gA = Vr — Vl appears only in R(D*). Vr and Vl receive constraints from 
both R(D) and R{D*). If new physics is established in both R(D) and R(D*) 
then the case of pure gA coupling is ruled out. The constraints on the complex 
couplings gy and qa are shown in the colored region of Fig. |6] (left) and (center). 
Fig. (0) (right) shows the constraints on the real couplings Vl and Vr. The real 
couplings are severely constraints by the recent R{D) and R(D*) measurements. 
The plot indicates that the data prefers either pure vector or pure axial vector 
couplings. Henceforth we will consider all three cases which include pure gy and gA 
complex couplings and real (Vl,Vr) couplings.. 
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Figure 6: The constraints on the complex coupling gy — Vr-\-Vl, 9a = Vr — Vl and 
the real (Vr, Vl) couplings at 95% C.L. 



3.1.2 Only gy coupling present 

In this section we consider only vector coupling gy = Vr + Vl- The NP vector 
coupling gy appears only in the amplitude A± and except for the forward-backward 
asymmetries, it does not significantly affect any other observables discussed earlier 
for the decay B — > D°*t~v t . One can see from Fig. [7] that the coupling g v can 




456789 10 456789 10 

? 2 [GeV 2 ] <? 2 [GeV 2 ] 



Figure 7: The left (right) panel shows the q 2 dependence of [Af B ]£>* ([v4^ b ]£>.) for 
the decay B — > D°*t~v t . The dashed lines show predictions for some representative 
values of gy. For example blue lines correspond to gy = 2.34e - * 3 ' 08 . 

enhance the magnitude of [A^bJd* up to 50% at low q 2 and it can have different 
zero-crossing point than the SM. In the SM, no zero-crossing is allowed for [A^]^*, 
however in the presence of gy, [A^ b ]e>* may have zero-crossing. Also, can 
reach up to 50% at low q 2 , and it can have either positive or negative sign. 
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3.1.3 Only g& coupling present 

In this section, we consider only pure axial vector coupling gy = Vr — Vl- In this 
case, except A± all other amplitudes depend on the new axial-vector coupling g& 
while the amplitude Ap is zero. Thus, the coupling g& does not significantly change 
the values for the polarization fractions of the D* meson and the r lepton from their 
SM predictions. In Fig. |8]we show the DBR and Rd* {q 2 ) in the presence of qa- The 
coupling gA can enhance DBR up to 0.4% GeV~ 2 at q 2 ~ 8.5GeV 2 , and Ro*(q 2 ) can 
be as high as about 0.9 at high q 2 . As shown in Fig. [9j the coupling gA can enhance 




? 2 [Gev 2 ] * 2 [Gev 2 ] 

Figure 8: The left (right) panel shows the q 2 dependence of DBR (Rpi*(q 2 )) for the 
decay B — > D°*t~v t . The dashed lines show predictions for some representative 
values of gA- For example, the red lines correspond to gA = 0.31e~* 2 ' 62 . 

the magnitude of [Afb]d* to about 50% at low q 2 and it can now have different 
zero-crossing point than the SM. In the SM, no zero-crossing is allowed for [A^Jd*, 
however in the presence of gA, [A'fbId* m ay have zero-crossing. Also, can 
take either positive or negative sign. The results are similar to the case where only 
gy coupling is present. 

3.1.4 Both Vl,r coupling are present and are real 

Finally, we consider the case where both Vl,r coupling are present and are real. In 
Fig. [TO] we show the DBR and Ro*{q 2 ) in the presence of both V L and Vr real cou- 
plings. These couplings can enhance the DBR upto 0.4% GeV -2 at q 2 ~ 8.5GeV 2 , 
and Rn*{q 2 ) can be increased to about 0.9 at high q 2 . 

One can see from Fig. [Til the couplings Vl and Vr can negatively enhance 
[^4fb]d* upto 50% at low q 2 and it can have a different zero-crossing than the SM. 
In the SM, no zero-crossing is allowed for [^4^b]d*, however in the presence of these 
new couplings [^4^ B ]z)* may have zero-crossing. Also, can reach upto 50% 

at low q 2 , and it can take either positive or negative values. 
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Only g A presents Only g A presents 




q 2 [GeV 2 ] q 2 [GeY 2 ] 

Figure 9: The left (right) panels of the figure show the q 2 dependence of [Aps^* 
([A-pbI-D*) for the decay B — > D°*t~v t . The dashed lines show predictions for 
some representative values of gA- For example the blue lines correspond to gA = 
2.34e- i3 - 08 . 




9 2 [GeV 2 ] ? 2 [GeV 2 ] 

Figure 10: The left (right) panel shows the q 2 dependence of DBR (Rr)*(q 2 )) for the 
decay B — > D°*t~v t . The dashed lines show predictions for some representative val- 
ues of (Vl, Vr). For example the red lines in correspond to (Vl, Vr) = (—0.97, 1.24). 



The polarization fractions of the D* meson are almost independent of the new 
couplings Vl and Vr. The tau lepton polarization fraction too does not depend on 
V L and V R . 

3.1.5 Pure Sl and Sr couplings present 

In this section we consider the scalar and pseudo-scalar couplings Sl^r. The combi- 
nation of the couplings Sr + Sl appears only in R(D), while Sr — Sl appears only 
in R(D*). If new physics is established in both R(D) and R(D*) then the cases of 
pure Sr, ± Sl couplings are ruled out. 

Hence, S# and Sl get constrained from both R{D) and R(D*). The constraints 
on the complex couplings Sr + Sl and Sr — Sl are shown in Fig. [T2] (left and 
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Only Vl.r Present Only V L , R present 
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Figure 11: The left (right) panel shows the q 2 dependence of [Afb].d* ([A^ s ]£>«) for 
the decay B — y D°*t~v t . The dashed lines show predictions for some representative 
values of gy. For example the blue lines correspond to (Vl, Vr) = (—1.02, 1.08) in 
the left panel and (Vl, Vr) = (—0.85, 1.21) in the right panel. 



middle). Fig. [T27 right) shows the constraints on the real couplings Sl and Sr. The 
real couplings are severely constrained by the recent R(D) and R(D*) measurements 
though the constraints are relatively weaker than those in the (Vl, Vr) case. To 
simplify our discussion we will take Sl,r to be real. 




Figure 12: The constraints on the couplings for complex (Sr + Sl) (left), complex 
(Sr — Sl) (center), and real Sl and Sr (right) at 95% C.L. . 



The combination of couplings Sr — Sl appears only in the amplitude Ap. In 
Fig. [13] we show the DBR and Rd* (q 2 ) in the presence of Sl and Sr couplings. 
These couplings can enhance the DBR up to 0.4% GeV -2 at q 2 ~ 7.5GeV 2 . Note 
that the peak of the DBR is shifted to low q 2 direction relative to the SM. The ratio, 
Rd*(q 2 ), can take the value of about 0.7 at q 2 ~ 7.5GeV 2 . 

The transverse forward-backward asymmetry [A^]d* is n °t sensitive to the 
Sl and Sr couplings . In Fig. [141 we show the effects of Sl and Sr on [Afb]d*, 
the polarization fractions (q 2 ) and P£ r (g 2 ). These couplings can positively or 
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Only S L , S present ^ S l>« P resent 




g 2 [GeV 2 ] q 2 [Ge\ 2 ] 

Figure 13: The left (right) panel shows the q 2 dependence of DBR (Rd*(q 2 ) ) for the 
decay B — > D°*t~j? t . The dashed lines show predictions for some representative val- 
ues of Sl and Sr. For example the red lines correspond to (Sl, Sr) = (1.11, —1.30). 



negatively enhance [A^b]^)* to about 30% at low q 2 and there can be different zero- 
crossing than the SM. Unlike the Vl and Vr case, the polarization fractions are 
sensitive to the Sl and Sr couplings. Due to the Sl and Sr couplings, F®* (q 2 ) 
can be as large as 0.85 at low q 2 , and it decreases to the SM value at high q 2 . The 
polarization fraction P L fT (q 2 ) can be as large as 0.5 and negative at very low q 2 and 
can have different zero-crossing than the SM. It increases to the SM value at high 
q 2 - 




456789 10 '456789 10 456789 10 



9 2 [GeV 2 ] ? ! [GeV 2 ] ? 2 [GeV 2 ] 

Figure 14: The left (center) and right panel shows the q 2 dependence of [A^s]^* 
(Ff *) and P£ T for the decay B -> D°* t u T . The dashed lines show predictions for 
some representative values of the new couplings Sl and Sr. For example the red lines 
corresponds to (Sl, Sr) = (—1.93, 1.73) in the left panel, (Sl, Sr) = (0.80, —0.9) in 
the middle panel and (Sl, Sr) = (0.61, —0.98) in the right panel. 
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3.2 B -> D°r _ zy T 

3.2.1 Only gy coupling present 

We now consider predictions for the various observables in B — y D°t~v t . The axial 
vector coupling does not contribute in this case and hence we consider only the 
coupling gy. Note that, the forward-backward asymmetry and the r polarization 
fraction are independent of the coupling gy. 

In the presence of gy , the q 2 dependence of DBR and Rn(q 2 ) for the decay B — y 
D t~v t are shown in Fig. [151 The DBR can increase upto 0.25% at q 2 ~ 7GeV 2 . 
The ratio Ro{q 2 ) is proportional to (1 + gy) 2 and increases with q 2 . 




<? 2 [GeV 2 ] ? 2 [GeV 2 ] 

Figure 15: The left (right) panel shows the q 2 dependence of DBR (R^(q 2 )) for 
the decay B —y D°t~v t . The dashed lines show predictions for some representative 
values of gy. For example the red lines correspond to gy = 0.68e* 1,27 in the left 
panel and gy = 0.62e iL20 in the right panel. 

Now we consider Vl,r real and independent. The q 2 dependence of DBR and 
Rr)(q 2 ) for the decay B —y D°t~u t are shown in Fig. [16J The predicted deviations 
from the SM are similar to the case with pure gy coupling. 

3.2.2 Pure Sl and Sr couplings present 

Finally, we consider the effect of Sl,r , taken to be real, on the observables in 
B — y D°t~v t . The allowed ranges for the real Sl and Sr couplings are shown as 
the colored region of Fig. IT2T right panel) . 

In Fig. [T7]we show the effect of Sl and Sr couplings on the DBR and Ro{q 2 ) in 
the decay B -> D q t-u t . The DBR increases up to 0.25%GeV" 2 at q 2 w 9.5GeV 2 . 
Note that the peak of the distribution in the DBR can be shifted towards high q 2 
relative to the SM. The deviation in the ratio Ro{q 2 ) increases with q 2 . 

Unlike the Vl,r case, the forward-backward asymmetry and r-polarization frac- 
tion in the decay B —y D°t~u t are very sensitive to the Sl and Sr couplings as 
shown in Fig. [JSJ In this case, can be either positive or negative, and 
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Only V L , R present ^ V l,R present 




q 2 [GeV 2 ] ? 2 [GeV 2 ] 

Figure 16: The left (right) panel shows the q 2 dependence of DBR (Ro{q 2 )) for 
the decay B — > D°t~v t . The dashed lines show predictions for some representa- 
tive values of new couplings Vl and Vr. For example the red lines correspond to 
(V L ,V R ) = (-1.09,-1.28). 



Only S L , R present ^ S l .* P resent 




9 2 [GeV 2 ] ? 2 [GeV 2 ] 

Figure 17: The left (right) panels of the figure show the q 2 dependence of DBR 
{Ro{q 2 )) for the decay B — > D°t~v t . The dashed lines show predictions for some 
representative values of new couplings Sl and Sr. For example the red lines corre- 
spond to (S L , S r ) = (0.51, -1.29). 

may have zero- crossing. In the SM , there is no zero-crossing for [Af_b]d • The 
T-polarization fraction can be negatively enhanced to more than 40% at low q 2 and 
may have the zero- crossing. 

4 Summary 

In summary, we have considered new physics explanation of the recent measure- 
ments in B — > D°*t~i? t and B — > Dt~v t decays by the BaBar collaboration. We 
considered an effective Lagrangian description of the new physics with four-fermi 
operators with vector/axial vector and scalar/pseudoscalar couplings. We consid- 
ered two cases, in the first case we considered only V/A couplings and in the second 
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Only S L , R present Only S L , R present 




4 6 8 10 4 6 8 10 

? 2 [GeV 2 ] « 2 [GeV 2 ] 



Figure 18: The left (right) panel shows the q 2 dependence of (P£ ) f° r 

the decay B — y D r~i? T . The dashed lines show predictions for some representa- 
tive values of new couplings Sl and Sr. For example the red lines correspond to 
(S L ,S R ) = (2.11,-0.04) in the left panel and (S L ,S R ) = (1.42,0.21) in the right 
panel. 

case we considered only S/P couplings. We found that the measurements in the 
decay B — > D°*t~v t strongly constrain the new physics couplings to be only pure 
vector or pure axial vector. Assuming a pure vector and pure axial vector com- 
plex couplings we calculated the differential branching ratios, the ratio Rf)*(q 2 ), the 
forward-backward asymmetries and the polarization fractions of the tau and the 
D* meson. We found that the pure complex vector couplings gy only affects the 
forward-backward asymmetries. The complex axial vector couplings, on the other 
hand, affects the DBR, Ro*{q 2 ) as well as the forward-backward asymmetries. The 
polarization fractions are not affected in presence of the gy,A couplings. When we 
considered real scalar couplings Sl,r , we found that all observables except the trans- 
verse forward-backward asymmetry for the D*, [Ap B ]D* were affected. Moving on to 
the decay B — y D°t~v t we found that the axial vector new physics coupling do not 
contribute to this decay. In the case of the vector coupling gy, the forward-backward 
asymmetry and the tau polarization fraction were not affected though the DBR and 
the ratio Rn{q 2 ) were affected. In the presence of the scalar couplings Sl,r, not only 
the DBR and the ratio Ruiq 2 ) were affected but the forward-backward asymmetry 
and the tau polarization fraction were also affected and were found to be very sen- 
sitive to the scalar couplings. The fact that different new physics couplings have 
different effects on the observables demonstrated that by measuring the various ob- 
servables it is possible to distinguish different models of new physics. Finally, if new 
physics is established in both R(D) and R(D*) then the cases of pure g^ = Vr — Vl 
and pure Sr — Sl couplings are ruled out as they contribute to only R(D*) and the 
pure Sr + Sl is ruled out as it contributes to only R(D) . 



17 



Acknowledgements 



We thank Gilad Perez, M. Papucci, M. Gonzalez-Alonso and J. F. Kamenik for 
useful discussion and comments. DG thanks Amol Dighe for discussion and encour- 
agement. This work was financially supported by the US-Egypt Joint Board on 
Scientific and Technological Co-operation award (Project ID: 1855) administered by 
the US Department of Agriculture and in part by the National Science Foundation 
under Grant No. NSF PHY-1068052. AD thanks the hospitality of the CERN theory 
group where the work was completed. 



Appendix 

A Kinematics 

The matrix element square for the decay B — > D^*Hui can be factorized into leptonic 
(L^) and hadronic (H^) tensors as 

\M{B -+ £>«^)| 2 = \(D^lu\c eff \B)\ 2 = L^H^. (22) 

The polar angle cos#; dependence of the leptonic and hadronic tensors L fll/ H fll/ 
can be evaluated using the completeness relation for the polarization four-vectors 
e(m = 0,±,t) [201 |2I]: 

J2 e»(m)? v {m')g mm , = <f , (23) 

m,m'=0,±,t 

where the tensor g mm i = diag(+, — , — , — ). The matrix element square reduces to 
\M(B -^D^lvi) | 2 = ]T L{m,n)H{m',n')g mm/ g nn/ , (24) 

m,m' ,n,n' 

where L(m,n) = L^ u e^{m)e* v {n) and H(m,n) = H^ u e*^{m)e v (n). The advantage of 
Eqj24]is that L(m,n) and H(m',n') are Lorentz invariant and so one can evaluate 
them in different Lorentz frames [21]. The leptonic tensor L[m,n] will be evaluated 
in the I — V\ center-of-mass (cm.) frame (q 2 rest frame), whereas the hadroic tensor 
H[m,n] in the B rest frame. 

In B rest frame, we choose the helicity basis e 

e(0) = -L(b D C.)|,0,0,- go ), e(±) = ±4=(0, ±l,-z,0), 

Vq 2 V2 

e(t) = _L( go ,0,0,|p D( .)|), (25) 

Vq 2 
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where q = (m B — m 2 D , + q 2 )/2m B and \p d m\ = \ l l 2 {m 2 B , m 2 D(it) , q 2 )/2m B . In this 
frame, the B and D^*' mesons four-momenta p B and p^w are 

p B = (m B ,0,0,0), p dW = {E dW ,0,0,\ PdW \), (26) 

where E D (*) = (m 2 B + m 2 D(t) — q 2 )/2m B . The momentum transfer q is q = p B —p D (*)- 
Further, one chooses the polarization vector of the D* meson as 

e(0) = —(\ PD .\,Q,0,E D .), e(±) = T ^(0,l,±i,0). (27) 

The leptonic tensor L[m, n] will be evaluated in the q 2 rest frame. In this frame, 
we choose the transverse components of helicity basis e to remain the same and 
other two components are taken as 

e(0) = (0,0,0,-1), e(t) = (1,0,0,0). (28) 

Let 9i be the angle between the D^*> meson and the r lepton three- momenta in 
the q 2 rest frame. We define the momenta of the lepton and anti-neutrino pair as 

pf = (E h pl sin9 h 0, -pi cos 9i) , 

Pv = (pi, -pi sin 9 h 0, pi cos 9i) , (29) 

where the lepton energy E\ = (q 2 + mf)/2y^ and magnitude of its three-momenta 
is pi = (q 2 — m 2 )/2y/q I . 



B B -> D*t~i7 t details 

The B — > D*t~i7 t differential decay rates for the lepton helicity A T = ±| are 
dT D *[\ T = -1/2] 



dq 2 d cos 9i 
dT D * [X T = 1/2] 
dq 2 d cos 9\ 



N\p D *\ 2\Aq\ sin 2 9 t + (\An\ z + \A ± \ 2 )(1 + cos 9^ - ARe[A\\Al] cos9 t 



m 



N\ PD .\-£ 2\A cos9 t - A tP \ 2 + (\A\\ \ 2 + \A±\ 2 ) sin 2 9 t 



The differential decay rate corresponding to the helicity A r = 1/2 vanishes for the 
light leptons (e, /x). 

The relevant form factors for the B — >■ D* matrix elements of the vector = 
cq^b and axial-vector A„ = 07^756 currents are defined as 



(D*(p D *,e*)\V, B{p B ) 
(D*(p D *,e*)\A^ B(p B ) 



*V P (J 

e p p D ,p B 



2iV{q 2 ) 

m B + mo* 

6* .q 

2m D *A (q 2 )—q fl + (m B + nio^A^q 2 ) 



(ms + mo* 



(Pfl + Pl?*)/ 



m 2 ^ — m 



(31) 



(30) 
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In addition, from Eq. ( I3TT) one can show that the B — >■ D* matrix element for 
the scalar current vanishes and for the pseudoscalar current reduces to 



(L>*(p D *,e*)|c 75 & B(p B 



2m D »A (q 2 



-e .q. 



(32) 



m b (n) + m c (/i) 

The expression of the hadronic helicity amplitudes for the B — > D*tu t decays 



arc 



An 



A, 



1 



[m 2 B - m 2 D * - q 2 ){m B + m D *)Ai( 



4m 2 B \p D *\ 2 
m B + m D * 



A 2 {q 2 ) (l-g A ) 



(m B + m^i(<? 2 )(l - 9 a) =F 2m ^ (g2) | Pd .|(i + ^) 



2m B |p D *|A (g 2 ) 
A = r^r (1 _ 9a) 



A f 



v 7 ? 
2m B \p D *\A (q 2 ) 

(m b (ii) + m c (fi)) 



9p 



C £> — > Dt v t results 

The B — > Dru T differential decay rates for the lepton helicity A T = ±| are 



(33) 



dT D [\ T = -1/2] 
dq 2 d cos 9i 
dT D [\ T = 1/2] 
dq 2 d cos 9i 



2N\p D \\H \ 2 sm 2 9 l , 

777^ 

2N\p D \^\H QC os9 l -H tS \ 2 . 



(34) 



The differential decay rate corresponding to the helicity A r = 1/2 vanishes for the 
light leptons (e, /i). 

The pseudoscalar form factors F + (q 2 ) and F (q 2 ) of the B —t D matrix elements 
are defined as 



{D{p D )\crfb B{p B 
(D(p D )\cb B(p B 



F + (q 2 ) Pb+ Pd . 



m B - m D ^ 



W) 



2n m B _-m^_ 



m B — m 



D 



m b (fi) - m c {jj) 
The helicity amplitudes are 



(35) 
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#0 



2m B \p 



F + (q 2 )(l + g v ), H t 



ni 



B 



m 



D 



F (q 2 )(l + g v ) 



m 2 B — m 2 D 
mb(fJL) - m c (/i) 



Fo(q 2 )gs. 



(36) 



D Form factors in the Heavy Quark Effective 
Theory 

In the heavy quark limit for the b, c quarks (m^ ^> Aqcd), both charm and the 
bottom quark in the hadronic current have to be replaced by static quarks h v ' )C and 
h Vt b, where v B = Pb / m B and v'^ t) = p^i*) l m D(*) are the four- velocity of the B and 
D(D*) mesons, respectively. The b — > c transition can be studied in the heavy quark 
effective theory (HQET). In this effective theory, the matrix elements of the vector 
and axial vector currents, and and , between bottom and charm mesons [23] 
are defined as 

(D(v')\V fl \B(v)) = ^b~^> [h+(w) + + h-(w) (v - 
(£>V,e')|l^|£(<;)) = i^m B m D * h v (w) e, ua ? e'*» v' a , 



(D*(v',e')\A ll \B(v)) = y/m B m D * h Al (w) (w + 1) e* - h A2 (w) e'*-vv 



-h Aa (w)e'*-v v' 



(37) 



where the kinematical variable w = vb-v D {*) = (m% + rn 2 D(t) — q 2 )/2mBm D (*). 

The form factors F + (q 2 ) and Fo(q 2 ) in Eq. fl35|) are related to the form factors 
h+(w) and h-(w) via 



R 



[l + w)R D 



Si(w) 



D 



(31 



where 



V x {w) 
Si(w) 



h + {w) 
h + (w) 



(1 + r) v 1 
(1 + r)(w- 1) 
(1 - r) (w + 1) 



h-(w) 



(39) 



Here Rp = 2^/mimD/ (mB+mc) and r = vtibIvhb- We will use the parametrization 
of the form factor V\{w) as given by [17] 



V^w) = 14(1) [1 - 8p 2 z + (51/?? - 10> 2 - (252p 2 - 84)z 3 



(40) 
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where z = (y/w + 1 — \/2)/(y/w + 1 + a/2). The numerical values of the free param- 
eters are [24] 

14(1)1^1 = (43.0 ± 1.9 ± 1.4) x 10" 3 , 

pj = 1.20 ±0.09 ±0.04. (41) 



For the form factor S±(w) we employ the parameterization as in [TT] . 

S^w) = [1.0036 - 0.0068(wl)±0.0017(wl) 2 ]Vi(w). (42) 
In the HQET, the helicity amplitudes in Eq. (1341) becomes 



H = m B (l±r) 



\ (1 +r 2 - 2rw) 



m B (l -r)y/r(w + 1) r . r , , m 2 B (l + r 2 - 2rw) i 

Hts = , ^^Si[w] [l + gv) - f — r~, r^\9s ■ (43) 

^(1 + r 2 _ 2rw ) 1 rni{m b {n) - mjji)) J 

Now we consider the form factors for the B — > D* matrix element in the HQET. 
The form factors h A .{w) are related to the form factors in Eq. (13~T]) [T0| \T7\ [25] in 
the following way, 

A 1 (q 2 ) = R D ,——h Al (w), A (q 2 ) = -^^h Al (w) , 

Mi 2 ) = ^-,H, v(f) = ^, M , (44) 

where Rd* = 2yfmBm* D / (mg + m£,). The w dependence of the form factors can be 
found in [TO], [T7] and the summary of the results are 

h Al (w) = /i Al (l)[l-8p 2 £±(53p 2 -15)^ 2 -(231p 2 -91)2 3 ] , 

R^w) = #i(l)-0.12(w- l)±0.05(u;-l) 2 , 

R 2 (w) = R 2 (l) + 0.11(w-l)- 0.06(w -l) 2 , 

Ro{w) = Rq(1) -0.11(w-1)+ 0.01(w - l) 2 , (45) 

where z = (y/w + 1 — \/2)/{\Jw + T+ v2). The numerical values of the free param- 
eters p 2 , 71^(1), -Ri(l) and i?2(l) are taken from[25j, 



/i Ai (l)|i/ c6 | = (34.6 ± 0.2 ± 1.0) x 10~ 3 , 

p 2 = 1.214 ±0.034 ±0.009, 

fli(l) = 1.401 ±0.034 ±0.018, 

i? 2 (l) = 0.864 ±0.024 ±0.008, (46) 
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and -Ro(l) = 1.14 is taken from [10]. In the numerical analysis, we allow 10% uncer- 
tainties in the i?o(l) value to account higher order corrections. 
In the HQET, the transversity amplitudes of Eq. ( 133]) become 



A, 



m B (l-r*)(w + l)-vAv, / \ 

// o ==— (™) 

y (1 + rj — 2r*w;) 



l)(l-J2a(«0) 



'1 - r* 



(1-<7a), 



= -m B y/2r*(w 2 - l)h Al (w)Ri(w)(l + g-y) , 
where r* = mn*/mB- 



(47) 
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